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^Pl Abstract 

P ■ 

\ Using QCD perturbation theory in NLO and light-cone QCD sum rules, we 

extract from the CLEO experimental data on the F 7 * 77r (Q 2 ) transition form factor 
constraints on the Gegenbauer coefficients ai and (Z4, as well as on the inverse moment 
{x^ 1 )^ of the pion distribution amplitude. We show that both the asymptotic and 
the Chernyak-Zhitnitsky pion distribution amplitudes are excluded at the 3a- and 
4<r-level, respectively, while the data confirms the end-point suppressed shape of the 
pion DA we previously obtained with QCD sum rules and nonlocal condensates. 
These findings are also supported by the data of the Fermilab E791 experiment on 
diffractive dijet production. 
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X 1 Introduction 

a: 

Perturbative QCD describes the short-distance interactions of quarks and gluons and 
can be applied to the description of hadronic reactions on account of factorization theo- 
rems. More precisely, one can calculate systematically perturbative kernels and associated 
anomalous dimensions that govern the evolution of hadron distribution amplitudes (DAs). 
These DAs parameterize hadronic matrix elements of quark-gluon currents and have to 
be determined by nonperturbative methods or extracted from experimental data. Re- 
cently, Schmedding and Yakovlev £Q have presented an analysis, based on light-cone QCD 
sum rules (LCSR) proposed earlier by Khodjamirian [2j and taking into account 0(a s )- 
corrections, of the high-precision CLEO experimental data [B] that allows to extract quite 
restrictive constraints on the first two Gegenbauer coefficients 02 and 04 which control the 
x-dependence of the pion distribution amplitude (7rDA). This sort of analysis was further 
extended and refined by us in ^ |3] with the aim to take more properly into account NLO 
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evolution effects of the 7rDA, to treat threshold effects of the effective strong coupling, and 
to estimate more carefully contributions resulting from (unknown) higher-twist effects. In 
addition, we derived directly from the CLEO data estimates for the inverse moment of 
the 7rDA, which is compatible with that obtained from an independent QCD sum rule, 
referring in both cases to the same low-momentum scale of the order of 1 GeV 2 . 

The results of our analysis, presented here, lead to the conclusion that the Chernyak- 
Zhitnitsky model jH| for the 7rDA in the plane (02, 04) is outside the 4cr-level, while the 
asymptotic DA, is excluded at the 3a level. In fact, the data seem to prefer end-point- 
suppressed DAs as those we have previously determined using QCD sum rules with non- 
local condensates [3]. These conclusions are further supported by contrasting the above 
mentioned 7rDAs with the E791 dijet data [7j following the convolution approach of Braun 
et al. [B]. Moreover, it was found [S] that the CLEO data are sensitive to the value of the 
average vacuum quark virtuality, limiting its value close to A 2 = 0.4 GeV 2 . 



2 What is the pion distribution amplitude (^(x, /i 2 )? 

The 7rDA is a central object in the deeper understanding of the pion microscopic structure 
in terms of quark and gluon degrees of freedom within QCD. This amplitude is defined by 
the matrix element of a nonlocal axial current on the light cone: 

(0\d(z) ltll5 E(z,0)u(0MP)) , = if*P* / ^e^rt// 2 ), (1) 



where gauge-invariance is ensured due to the Fock-Schwinger string E(z, 0) = Ve %9 ^ A v.( T ) dT>1 
and ip^ w ' 2 (x, /i 2 ) is symmetric with respect to x <-» x (x = 1 — x) and is normalized to 
unity, whereas fi 2 denotes the normalization scale. Fig. ^ visualizes the light-cone struc- 
ture of the 7rDA. There are also 6 pion DAs at twist-4 level, four of them contributing to 




Figure 1: <p w (x; jj 2 ) - light-cone amplitude for the transition n — > u + d. 



the 7*7 — > 7r-transition as a twist-4 correction, whose value is parameterized by the scale 
~ 0-19 GeV 2 . In what follows we will speak mainly of the twist-2 7rDA and for the 
sake of brevity we will omit its superscript Tw " 2 referring to it simply as (p n (x; /1 2 ). 

Due to vector current conservation, the solution of the ERBL evolution equation [5J EH] 
(in LO approximation) in the asymptotic limit is <p n (x; fi 2 — > 00) = ip As (x) = 6x(l — x). A 
particularly convenient way to represent the 7rDA is to use its 1-loop eigenfunctions, viz., 
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the Gegenbauer polynomials 



^3/2 



^(s; /i 2 ) = ^ As (*) 1 + a 2 (n 2 )Cr(0 + a^Cf + ... (2) 

t;=2x — 1 



with C n (0 being the Gegenbauer polynomials and the ellipsis denoting still higher-order 
eigenfunctions than displayed. In this representation all the dependence of (p n (x; fi 2 ) on ji 2 
is concentrated in the coefficients a n (/z 2 ) due to the fact that the 1-loop evolution kernel has 
a factorized structure Vi_i oop (x, x'\ a s ) = [a s / (47i)]V (x, x'). In the NLO approximation the 
eigenfunctions of the evolution kernel inevitably depend on a s and therefore on a 2 . Note 
that because of the symmetry in x «-> x, only even Gegenbauer polynomials contribute. 

The high precision of the CLEO data provides the possibility to extract these important 
theoretical parameters (a 2 and 04) directly from experiment. But before we turn to this 
subject, let us first give some brief exposition of the theoretical method to determine the 
7rDA within the QCD sum-rule approach. 



3 QCD sum rules with nonlocal condensates 

To model the nonlocality of the QCD vacuum, we assume (q(0)q(z)) = (g(0)g(0))e~' z2 ' A «/ 8 , 
and similar expressions for other nonlocal condensates (NLCs), where a single scale param- 
eter A 2 . = (k 2 ) was introduced in order to characterize the average momentum of quarks 
in the QCD vacuum jTTj: 

r 0.4 ± 0.1 GeV 2 from QCD SRs [12] ; 
A 2 , = 1 0.5 ± 0.05 GeV 2 from QCD SRs [T3|; 

( « 0.4- 0.5 GeV 2 from lattice QCD OUS]. 

The correlation length A^ 1 ~ 0.3 Fm ~ p-meson size represents the width of the NLC 
at small distances. Let us mention that for very large distances {z ^> 1 Fm 15J one 
may assume another form of the condensate, given by jTH] (q(0)q(z)) ~ (<?(0)g(0))e~' z ' A 
at \z\ ^> 1 Fm (with A ~ 450 MeV). This behavior is of no importance in the problem 
under investigation. 

In [4J we have determined all coefficients up to order n — 10 using QCD sum rules 
with nonlocal condensates. It turned out that all coefficients beyond n = 4 are very small 
so that for practical purposes it suffices to model the 7rDA using only a 2 and a 4 . So, the 
NLC QCD sum rules produces a whole bunch of self-consistent 2-parameter model DAs 
(see Fig. 2a.) at /1 2 ~ 1 GeV 2 : 

<p n (x) = <p as (x) [l + a 2 C 3 2 /2 (2x - 1) + a 4 C 3 J 2 (2x - 1)] (3) 

with the best-fit model (bold-faced on Fig. 2a) defined by the parameters 

aY' = +0.188; < f - = -0.130. (4) 

The admissible regions for the parameters a 2 , a 4 of the 7rDA are presented in Fig. 2b as 
shaded slanted rectangles and are shown for different values of A 2 . Fig. 2a demonstrates 
the most striking feature of our type of 7rDAs: their end-points (i. e., x — ► and x — ► 1) 
are strongly suppressed, the suppression being controlled by the quark vacuum virtuality 
A 2 . Both the asymptotic and the CZ 7rDAs are not end-point suppressed, as we have 
quantitatively shown in |4,. Our models demonstrate by a precedent that the common 
statement "two-humped tiDAs are end-point concentrated" is wrong. 
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4 7*7 ~~ 71-1 Why Light-Cone Sum Rules? 

For Q 2 ^> m 2 , q 2 <C m 2 , pQCD factorization does not help because it is valid only 
for leading twist and therefore higher twists become important ^7j. The reason for this 




Figure 3: Left part demonstrates the regime when pQCD description is valid; right part 
makes explicit why LCSR should be applied. 

failure can be understood by recalling that if q 2 — > 0, one needs to take into account the 
interaction of a real photon at long distances of the order of 0(1/ >/q 2 ~), as the following 
Fig. 3 illustrates. To account for long-distance effects in a perturbative QCD treatment, 
one needs to introduce the light-cone DA of a real photon. 

To this end, Khodjamirian [2] has shown that light-cone QCD sum rules (LCSR) 
effectively account for the long-distance effects of a real photon by using quark-hadron 
duality in the vector channel and an appropriate dispersion relation in q 2 \ namely, 

F ^ iQ > q) = nJ ^TT 2 exp [h^J ds+ n J so ~JT^f ' (5) 

where so — 1.5 GeV 2 is an effective threshold in the vector channel and the Borel parameter 
M 2 takes values in the range 0.5 — 0.9 GeV 2 . Then, the real photon limit (q 2 — > 0) becomes 
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safely accessible. Here p(Q 2 , s) = lmF^, n (Q 2 , —s) includes contributions from both the 
leading twist 7rDA as well as the twist-4 one. The latter is characterized by the twist-4 
scale parameter S^ w _ 4 . This theoretical ground was extended by Schmedding&Yakovlev 
(SY) to the NLO accuracy [Tj. 

5 Results from nonlocal QCD sum rules vs CLEO 
constraints 

In jSj we improved the SY analysis based on LCSR (jHJ) by taking into account ERBL NLO 
evolution for the 7rDA and the exact NLO running of a s (Q 2 ). The established relation 
~ has been also involved in the analysis. As Fig. 4a shows, we obtained 

reasonable agreement with established in this approach constraints just for the value of 
X 2 q = 0.4 GeV 2 . 




Figure 4: Comparison of theoretical predictions of nonlocal QCD sum rules for X 2 = 
0.4 GeV 2 and the CLEO data constraints obtained in the LCSR approach. Left (a): 
Previous [5] results. Right (b): New BMS [IS] constraints. Here: ♦ = asymptotic DA, 
X = BMS model, ■ = CZ DA, + = best-fit point, [H] and ♦ 2I| = instanton models, 
T = transverse lattice result [21]. All values are evaluated at /ig Y = (2.4 GeV) 2 . 

More recently we have refined this extensive analysis in several respects, notably, 
by obtaining from the CLEO data direct estimate for the inverse moment of the 7rDA that 
plays a crucial role in pion electromagnetic/transition form factors and by verifying the 
reliability of the main results of the CLEO data analysis quantitatively. We also refined 
our error analysis by taking into account the total uncertainty of the twist-4 contribution 
and treated the threshold effects in the strong running coupling more accurately. The 
main upshot of this investigation is presented graphically in Fig. 4b, where the values 
A 2 = 0.4 GeV 2 and <5f w _ 4 = 0.19(4) GeV 2 have been employed. One can see that even 
with a 20% uncertainty in the twist-4 contribution, the CZ distribution amplitude (■) is 
excluded - at least - at the 4cr-level, while other well-known models (*&, ♦ and T) with 
shapes more or less close to the asymptotic one (♦) are excluded at the 2cx-level. 

These findings are further supported by extracting the inverse moment of the 7rDA 
from the CLEO data in a two-Gegenbauer model, a 2 + a 4 = (x _1 )^. xp (/iQ)/3 — 1, at the 
low scale /i 2 , = 1 GeV 2 . The obtained constraints are presented in Fig. 5b. One should 
compare them with the theoretical model-independent estimate of the inverse moment 
(x _1 )^ R (/ig « 1 GeV 2 ) = 3.28 ± 0.31, obtained in the special NLC QCD sum rule using 
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Figure 5: Left (a): The inverse moment (x~ 1 )^ shown as a function of the Borel param- 
eter M 2 from the special NLC QCD sum rule at the scale [3]; the light solid line is the 
estimate for (x _1 )^ R ; the dashed lines correspond to its error-bars. Right (b): The result 
of the CLEO data processing for the quantity (x _1 )° xp /3 — 1 at the scale /i 2 , ~ 1 GeV 2 in 
comparison with three theoretical models from QCD sum rules, CZ, BMS, and (a). The 
thick solid-line contour corresponds to the union of 2cx-contours, while the thin dashed- 
line contour denotes the union of lex-contours. The light solid line with the hatched band 
indicates the mean value of (x _1 )^ R /3 - 1 and its 

error bars in part (a). 



again A 2 = 0.4 GeV 2 [22 E], see Fig. 5a. Noteworthily, these constraints match each other 
and both of them comply with the value ^(a;^ 1 )^ — 1 = 0.24 ± 0.16 found in fS] from 
a LCSR analysis of electromagnetic pion form factor. From Fig. 5b it is evident that 
again both the asymptotic 7rDA and the CZ model are far outside the region of the CLEO 
experimental data. 



6 E791: Diffract ive dijet production 

The Fermilab group E791 proposed [7j to exploit experimentally the ideas on dijet diffrac- 
tive dissociation suggested in [23] and further developed in [231221 IE] Braun et al. [E] 
have used a convolution-type approach to account for hard-gluon exchanges, represented 
diagrammatically in the left part of Fig. 6. Following this convolution procedure (hav- 
ing also recourse to [2I|), and ignoring the distortion of our predictions caused by the 
detector acceptance, we found the results displayed in the right part of Fig. 6, making 
evident that, though the data from E791 are not that sensitive as to exclude other shapes 
for the pion DA (asymptotic and CZ model), also displayed for the sake of comparison, 
they are relatively in good agreement with our predictions. Especially, in the middle x 
region, where our 7rDAs - the shaded strip - has the largest uncertainties (see Fig. 2a), 
the predictions are not in conflict with the data. However, before this data set can be 
used for a quantitative comparison, its inherent uncertainties have to be removed. 

It is again worth emphasizing that because our model distribution amplitudes - exem- 
plified by the BMS model - are end-point suppressed (see Fig. 7), they are not affected 
by the poor accuracy of the E791 experimental data in these regions. 
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Figure 6: Left: Diffractive dijet 7Tv4-production in the E791 experiment with q\ ~ 4 GeV 2 
and s ~ 1000 GeV 2 . Right: Asymptotic DA (solid line), CZ DA (dashed line) and BMS 
bunch (shaded strip) in comparison with E791 data. Corresponding x 2 are: 12.56, 14.15 
and 10.96 (the last for BMS model with A 2 , = 0.4 GeV 2 ). 

7 Conclusions 

Thanks to the recent high-precision CLEO experimental data [3], we can answer more 
questions of nonperturbative QCD than a couple of years before. On the theoretical side, 
the method of QCD sum rules with nonlocal condensates [TTJ |22J UIH H| provided a tool 
to determine more precisely than before a bunch of candidate DAs for the pion that are 
endpoint-suppressed due to a rather large QCD vacuum quark virtuality A 2 . On the other 
hand, the method of light-cone sum rules [3 E E] enables us to access the pion-photon 
transition form factor when one photon becomes real. Taking these theoretical approaches 
in conjunction, we were able to analyze the CLEO data at the NLO level in order to derive 
restrictive constraints on the first two Gegenbauer coefficients <22 and 0,4, which control 
the x-dependence of the 7rDA. These parameters allow the reconstruction of the 7rDA and 
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Figure 7: Comparison of different DA 
curves aiming to illustrate the end-point 
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can be further tested against other experimental data, like those collected in the dijet 
production Fermilab experiment E791. The main conclusion is that both the CZ model 
as well as the asymptotic 7rDA are excluded — at least at the 2a level — by the CLEO data, 
while the two-humped end-point suppressed BMS distribution amplitude with a value of 
\ 2 q ~ 0.4 GeV 2 is in a good agreement with the CLEO data and not in contradiction with 
the E791 data. 
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